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A formulation of variational principles in terms of functional integrals is proposed for any type 
of local plastic potentials. The minimization problem is reduced to the computation of a path 
integral. This integral can be used as a starting point for different approximations. As a first 
application, it is shown how to compute to second-order the weak-disorder perturbative expansion 
of the effective potentials in random composite. The three-dimensional results of Suquet and Ponte- 
Castafieda (1993) for the plastic dissipation potential with uniform applied tractions are retrieved 
and extended to any space dimension, taking correlations into account. In addition, the viscoplastic 
potential is also computed for uniform strain rates. 
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I. INTRODUCTION 

In the last decade, studies of disordered non-linear composites have forced new homogeneization methods in Me- 
chanics. Most such methods rely on rigorous variational principles, from which exact optimal bounds can be deduced 
(Ponte-Castaheda and Suquet, 1998). These principles are either standard minimum energy principles, or more refined 
approaches using a linear reference material, be it homogeneous (Hashin and Shtrikman, 1962ab; Talbot and Willis, 
1985), or even heterogeneous (Ponte-Castaheda, 1991; Suquet, 1993). 

In the study of non-linear dielectric media as well, few exact results were previously available, save for the case 
where the non-linearity in the response was treated as a perturbation of the linear behavior (Stroud and Hui, 1988). 
In a recent pioneering work, Blumenfeld and Bergman (1989) computed the effective dielectric permittivity of a 
strongly non-linear random medium to second order in a weak-disorder expansion; cf. also Bergman and Lee (1998). 
The perturbative expansion was carried out by use of the Green function associated to the linear problem. An 
analogous approach was subsequently undertaken by Suquet and Ponte-Castaheda (1993) on strongly non-linear 
elastic composites. These results were then further extended to composites with inclusions of complex shapes (see, 
eg., Ponte-Castaheda and Suquet (1998) for a review). 

Perturbative results are important as testing-benches for self-consistent effective-medium theories. Linear self- 
consistent formulae (Budianski, 1965) are known to give qualitatively correct predictions even for a high contrast 
between the constituents of the random medium. They account for, e.g. the existence of a "rigidity threshold" , 
i.e. a volumetric concentration of voids above which a porous linear elastic material loses its rigidity (Sahimi, 1998). 
In addition to the fulfillment of exact non-linear bounds (Gilormini, 1995), criteria of acceptability for a non- linear 
self-consistent formula include the recovery of second-order perturbative results, and also constraints on the critical 
behavior near the rigidity threshold (Levy and Bergman, 1994) where field fluctuations are enhanced. In Mechanics, a 
recent self-consistent effective-medium theory (Nebozhyn and Ponte-Castaheda, 1997) reproduces exact perturbative 
results to second order, but its critical fluctuations have not been studied yet. An equivalent theory has been proposed 
in electrostatics (Ponte-Castaheda and Kailasam, 1997). These result were derived from a variational principle. 
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As a complement to classical means, functional methods originating from field theory (for a book on functional 
methods see, e.g., Kleinert, 1995) and from the statistical physics of disordered systems (Mezard et al., 1987) were 
recently harnessed to contribute to effective-medium studies (Barthlmy and Orland, 1993; Parcollet et al., 1996). They 
have not yet spread in the mechanics community. Their major interest is that they not only provide a convenient 
workhorse for implementing variational principles under various constraints, but that they also lend themselves to 
approximations resulting in self-consistent effective-medium formulae (Budianski, 1965; Parcollet et al., 1996), some 
of them being of a new type (Pellegrini and Barthelemy, 1999). With such methods, the problem of the minimization 
of the potential becomes equivalent to the computation of a functional integral, as shown below. This single integral, 
which may be written down for all types of local potentials, can then be approximated using standard tools (saddle 
point methods, perturbative expansions, etc.). 

The purpose of this article is to introduce in detail the specific mathematical apparatus needed to apply functional 
methods to variational principles in mechanics. As a first application, we re-derive the result of Suquet and Ponte- 
Castaheda (1993) for the weak-disorder expansion of a strongly non-linear viscoplastic composite assumed to be of 
the Norton type. We extend it to any space dimension, and show how to take correlations into account. 

We adopt an approach (Barthelemy and Orland, 1998) which utilizes the so-called "replica method" (Edwards and 
Anderson, 1975) developed in the framework of spin glass theory (Mezard et al, 1987). The replica method enables one 
to average over the disorder in a non-perturbative way, but is by no means compulsory in order to exploit the functional 
formulation. It is however a most natural route towards extensions to self-consistent formulae. Specific applications 
to self-consistent formulae for linear or non-linear materials (including the EMT) will be presented elsewhere. 

The paper is self-contained, and is organized as follows. Section || introduces the variational formulations of the 
problem and explains how to cast them under a functional form, either starting from the local viscoplastic potential 
(expressed in terms of the stress tensor) or from the plastic dissipation potential (expressed in terms of the strain rate 
tensor) . We show that minimizing a potential amounts to computing the statistical average of the logarithm of some 
suitably defined partition function. The computation is carried out with the replica method, explained in t he tex t. 
In Sec. [II, we specialize to the perturbative expansio n, begin ning with the plastic dissipation potential (Sec. Ill A 1 ) . 



The results are then applied to the Norton law (Sec. Ill A 2 ). We retrieve the result of Suquet and Ponte-Castaneda 
(1993), generalized to any space d imen sion. The perturbative calculation of the viscoplastic potential, slightly more 



involved, is presented next, in Sec. Ill B . In both cases, the disorder is assumed to be site- disorder, the local potentials 
are statistically uncorrelated from point to point, for simplification purposes. This restriction is overcome in section 
IV, where spatial correlations (not necessarily isotropic) are introduced. Before concluding in Sec. VI, we comment 
on the range of applicability of the method (Sec. |v|). 

Technical details are left to appendixes. Appendix ^ is a brief reminder of the use of fourth-rank tensors in 
mechanics, and defines various fundamental fourth-rank tensors employed hereafter. Various algebras and sub-algebras 
encountered in the paper are also defined in this appendix, which ought be read before the reader goes through our 
calculations. The determinant (det), trace (tr) and inversion ( _1 ) operators will always be indexed by a label indicating 
in which algebra or sub-algebra of operators they act. Useful formulae for restricted gaussian integrals over vectors or 
second-rank tensors, which we could not find in the literature, are derived in Appendix |b| Inversion and determinant 
formulae for matrices in the replica space are given in Appendix Evaluations of functional integrals are relegated 
to Appendix The Legendre transform between the perturbative expansions of the viscoplastic and of the plastic 
dissipation potential is examined in Appendix [e]. 

Before going on, we state our notational conventions: unless otherwise indicated, scalars are denoted by regular 
typefaces (eg. a or A); vectors are denoted by bold typefaces (a); tensors of rank two by sans-serif typefaces (a) 
[the exceptions are <r and £ which appear in boldface], and tensors of rank four by capital Blackboard typefaces 
(A). Tensor indices are denoted by roman letters (e.g. A^), whereas replica indices are Greek letters. The Einstein 
summation convention on repeated indices is used, except when the indices are underlined (An = tr A, whereas An 
denotes the i th diagonal element of A). Dyads (i.e. tensor products of vectors) are denoted omitting the tensor product 
operator: A = ab is the dyad with components Aij = aibj. 



II. THE PROBLEM AND ITS FORMULATION WITH REPLICAS 



A. Variational principles and effective potentials 



We adopt here the presentation of Ponte-Castaheda and Suquet (1998). Let <r(x) be the local stress (tensor) field, 
and d(x) be the local strain rate (tensor) field deriving from the local velocity (vector) field v(x). The material 
domain, of volume V, is denoted by 0. Then the equilibrium equations on er read, for x e il (the upper index * 
denotes the transpose): 
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V.er = 0, (2.1) 

cr = t cr. (2.2) 

In addition, 

d = ^[Vv + t (Vv)]. (2.3) 

The local constitutive relations are expressed by means of either the viscoplastic potential ip x or the plastic dissipation 
potential </> x , as: 

- = ^(d), d = — (<r). (2.4) 

The subscript x indicates that the potentials may vary from point to point in the material, according to the local 
material properties. Both local potentials are usually convex functions of the fields, and are linked by the (Legendre) 
duality relation 

</> x (d) = max [tr :d — ip x (<r)] . (2.5) 

(T 

Classically, two types of boundary conditions are considered: 

1) for boundary conditions of uniform traction (ut) on the boundary <9f2 of the material, the heterogeneous material 
is shown to be macroscopically described by a pair of effective potentials \? and $ such that 



* ut (S) = min Vx(t), (2.6a) 

°7\ o-.n=£.r, on 



. an 



$ nt (D) = min x (d). (2.6b) 

v/d=D 

(the over line denotes a volume average on f2, and n is the outward normal to f2). 

2) for boundary conditions of uniform strain rate (us) D on the boundary, the effective potentials are: 



# U8 (£) = min ^ x (<r), (2.7a) 



$ US (D) = min x (d). (2.7b) 

v/v— D.x on 

For both types of boundary conditions, the effective potentials $ and \& are shown (Suquet, 1987; Willis, 1989) to 
be (Legendre) dual functions such that 

• (E) + .(D) -3:0, E . (2.8) 

In both cases X = <r and D = d. Depending on the boundary conditions considered, the latter equalities are either 



definitions, or theorems (Hill, 1963). Therefore only one equation in each pair (2.6) or ( |2.7| ) is sufficient to characterize 
the effective homogeneized material. 

The pro blem considered in this article consists in computing the effective potentials for disordered materials. Eqs. 



(2.6), (2.7) refer to one sample (one particular realization of the disorder). A key assumption is that the effective 
potentials are self-averaging, i.e. that they do not depend on the sample at hand in the so-called "thermodynamic 
limit" , where the size of the sample goes to infinity. In other words, the material contains in the various regions of its 
bulk all the possible realizations of the disorder. Therefore, rather than carrying out the calculation of the effective 
potentials for one particular configuration of the disorder (which is impossible), we arrive at the correct result by 



considering statistical averages. For instance, in the limit V — > oo we can as well write, instead of (2.6a): 



tt ut (E) = (*ut(E)> = ( min ■ (2-9) 

w{ ,r.r >:'„"'■', I 

Note that the statistical averaging operator (.) and the infimum operator do not commute (inverting the operators 
would lead to the incorrect result that the effective potential is the average of the local potentials). 
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B. Functional representations 



The minimization problem is reformulated by extending to functionals the following straightforward result for scalar 
functions: 



[a,h]M = - lim Jlog f dxe-W'l 

P J a 



(2.10) 



Note that this formula holds even if the minimum is met at several points in the interval. 

Each infimum ( |2.6| ) or ( |2.7| ) is obtained for the state of the fields <r or d which minimize the potentials, among 
all the possible states fulfilling the equilibrium and boundary constraints. In statistical mechanics, this problem is 
analogous to that of finding the lowest energy state Eq of some hamiltonian H[s] functionally depending on some field 
s(x). In this case, Eq is given by 



E 



lim - log Z, 



Z = tr e 



-0H 



(2.11a) 
(2.11b) 



where the trace operator tr here denotes the sum over all the allowed states of the system, i.e. the configurations of 
s. Indeed, the partition function Z is obviously dominated by exp(— /3Eq) in the limit /3 — > oo, where Eq = H[sq] and 
so is the configuration of s that minimizes H (/? is the reciprocal of the temperature in statistical physics). 

The means of summing over all the possible states of a continuous field subjected to constraints is to use functional 
integrals (also termed path integr als) w ith sui tably defined measures. With such a formalism, it is more convenient to 
use the variational formulations ( |2.6b ) and ( [2.7a ) where the boundary conditions are implemented through volume 
averages. We therefore write (with the additional statistica l ave rages): 

1) for uniform tractions [d is computed in terms of v by J2.3|)1, 



*ut(D) = - 

where S is the Dirac distribution and 
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lim j- (log Z ut ) , Z ut = \Vw S(d - D) e-^W 

(3— too fj 



xGfi 

S(d-D) = Y[S(d lJ -D l3 ); 



For brevity, the constrained measure in (2.12) will be denoted by 

£> ut v = 2Va(d-D). 

2) for uniform strain rates, 

1 



*„.(£) = - hm -(logZ u 

p—foc p 



with 



xGfi 

d s a( X ) = 2 d ( d - 1 )/ 4 [[|*(^(x) - ^(x))][n^W], 

Kj 



(2.12) 

(2.13a) 
(2.13b) 

(2.14) 

(2.15) 

(2.16a) 
(2.16b) 
(2.16c) 
(2.16d) 



5(W-S) = ]JS(a ij -'E ij ). 

i<j 

The integration measure d s cr(x) allows one to integrate over all symmetric stress tensors (cf. Appendix |b|). The 
normalization factor 2 d ( d_1 )/ 4 in d s er(x), where d is the space dimension, is of no importance but has been included 
for consistency with Appendix pi We set: 



X> us er = V s (t5(V.(t)5(W - S). 



(2.17) 
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C. The replica method 



Hereafter, we focus on Z ut , but the following considerations apply mutatis mutandis to Z us . The statistical average 
of log Z ut is a complicated quantity to evaluate directly. The replica method consists in writing 



(logout) = lim -((^ t ) - 1), 
r— >0 r 



(2.18) 



and in computing (Z^ t ) with r an integer (in the replica literature, the number of replicas, r, is usually denoted by n. 
We do not use the latter in order to avoid confusions with the Norton exponent). At the end of the calculation, we 
take the limit r — > 0, assuming the analytical continuation is legitimate, which we cannot prove in general. Therefore 
introducing r replicas v Q of the field v, the partition function 



^ t = ri (/ 

a=l KJ 



3 -/30 x (dT) _ 



n^tv 



(2.19) 



is that of a system made of r replicas of the initial system. The greek labels (which exclude /3, this symbol being 
devoted to the reciprocal of the "temperature" ) are replica labels running from 1 to r. The statistical average therefore 
is 



m = 



J n ^ 



-/3E; = i0=c(dT) 



(2.20) 



For simplification purposes, we shall limit ourselves to site-disordered materials where the local potentials are 
statistically uncorrelated from point to point (we explain in Sec. IV how to overcome this restriction). By definition 
of site disorder, we have for any function F the property that 



(F (Vx(<t(x))) F (^ y (<x(y)))) = (F (Vx(<x(x)))) (F (Vy(«r(y)))) if x ± y. 



(2.21) 



For this type of disorder, the average in ( |2.20| ) is readily carried out. Replacing the volume integral by a Riemann 
sum J dx — > i>X} x s w here v denotes the unit "infinitesimal" volume element of the theory (here the size of v also is 
the statistical correlation radius within which the potential is constant), one obtains 



Q = l 

a=l \ : 

ri ^ a n 



a=l 



-(v0/v)T,: =1 4>^''y 



This finally yields: 



J ni? ut v«exp [/^log 

a—1 



'(vfJ/V)Z:, =1 Md-') 



(2.22) 



(2.23) 



where we reintroduce the integral notation in the exponent. This is the usual shorthand notation used in functional 
integral methods and has to be un ders tood a s the limit of a Riemann sum. 

This e xpress ion, completed by ( 2 . 1 g ) and ( 2.12 ) , constitutes an exact functional representation of the variational 
problem (2.6t) for the particular type of disorder under study, and is valid for any form of x (d). In other words, 
we recast a minimization problem under a computational form. Computing this integral exactly is equivalent to 
solving the minimization problem. We note that before averaging, the replicas are independent from one another, 
but are coupled after the average is carried out. Here appears an important feature of the replica method: the use 
of replicated fields allows one to average over disorder by transforming the disordered problem into a non-disordered 
one, but with couplings (if there were no couplings, averaging the logarithm would be a trivial matter). We are now 
left with a complicated field-theory problem without disorder, consisting of r coupled fields, and we have to resort to 
approximations. In the next section, we present a perturbative calculation of this integral. 
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An important simplification arises from the fact that (Z* t ) always shows up in the form (Z^ t ) = (Cf3 c ) r f(r, (3), 
where /(0, x) = 1 and is a non-homogeneous function of x, C is a /3-independent constant, and c is an exponent linked 
to d, the space dimension. Therefore 



log(C/3 c ) + i log f(r,f3) + 0(r) 



(2.24) 



and the multiplicative factor (Cf3 c ) r does not contribute to the limit (3 — > oo. Hereafter, such multiplicative factors 
will therefore systematically be dropped out of the calculations. This will be indicated by the proportionality symbol 
oc appearing instead of the equal sign in equations. 

III. PERTURB ATIVE CALCULATION 

We detail in this section the perturbative calculation for $ first, since it is simpler than that for The latter 
is examined next, once the reader has become acquainted with the method. In the following, we assume that the 
c?-dimcnsional material is rigid- viscoplastic, described by the dual [in the sense of (2.5)] Norton potentials (1 < n < oo, 
m = l/n). 



x (d) =6 m d™+ 1 /{m + l) 
tr 2 (d) = 



^ x (<T)= Wnf T" +1 /(n+l), 



with 



£ /<7™. 



The equivalent Mises norms cr eq and d e „ are: 



d- 1 



d 2 



W 2 ) 



1/2 



1 2 



1/2 



(3.1) 



(3.2) 



(3.3) 



These are the d-dimensional counterparts of the usual three-dimensional definitions (for d — 3, we recover the factors 
2/3 and 3/2). The deviatoric part of the stress or strain rate tensor has been denoted by a prime: a' = a — tr 2 (a)l/d, 
where a = er or d. 

The additional incompressibility constraint tr 2 (d) = V.v = 0, which characterize incompressible materials, can as 
well be accounted for by adding a term (K/2) tr 2 (d) 2 to 4> x (d), and by letting the constant K — > oo at the end of the 
calculation. Instead, we shall implement it directly by multiplying the measure T> ut v by a factor d(V.v). Both ways 
actually amount to the same. 

The constitutive parameters gq and e'o decribe the material and are random functions of the position variable x, 
uncorrelated from point to point, according to our site-disorder hypothesis. The exponent n is taken to be constant 
in the material, for simplification purposes (but this is not a limitation of the method) . 



A. The plastic dissipation potential 



1. General framework 



Integral Q2.23D is an exact expression for the effective potential. Its perturbation expansion is considered hereafter. 
The measure V ut v a contains the term 5(6 — D). This means that the main contribution to this integral is coming 
from fields around D. Indeed, if we suppose that only the constant field d = D (i.e. v = D.x such that tr 2 (D) = 0) 
contributes to the integral, we find the trivial lowest order result 



3„t(D) = (0 X (D)). 



(3.4) 



We shall now seek the first corrective term, and compute ( [2.23 ) by integrating over fields v such that d differs little 
from D. The perturbation expansion is obtained with Laplace's method for asymptotic expansions of integrals (Bender 
and Orszag, 1978). We therefore write 



d = D + e 



(3.5) 
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and expand the logarithm in ( 2.23] ) in powers of e. The linear and quadratic terms are then kept inside the main 
exponential, and the rest is further expanded in e. We start from the expansion: 



</> x (D + e) = X (D) + WLea + X Meji e lk + C7(e 3 ), 



with 



<-(x) 



ddij (x) 



d 2 ^ 



d=D 



ddij(x)ddki(x) 



(3.6) 



(3.7) 



d = D 



With a change of integration variables such that v now denotes the correction to D.x, the expansion of (Z^ t ) takes 
the form (taking into account the incompressibility constraint): 



x exp 



« / i[v ut v a 5(V.v a ) [l + 0(e 3 )] 
-/3(W') e : J2 ^ - \ ^ E / dx e ° W :M ° 7 : e7 W 



(3.8) 



The functional measure now is 2? u tv = XV<5(e), and e = [Vv + (Vv)]/2. Besides, we have introduced the matrix 
MP 7 = Mi^ Q7 - (3(v/V)M 2 , with 



h = <W")e 

h = (W'W')e - (W')e(W') 



The notation (-) e stands for the weighted average 



(3.9a) 
(3.9b) 



(3.10) 



Only the leading term in r is needed in ( [3.10 ) [a handwaving argument for this is that the sums over replicas in ( |3.8| ) 
are at least proportional to r: the contribution is null for r = 0. They therefore already provide the linear terms in r 
required for evaluating ( 2.18| )], and we shall therefore in practice identify (.) e with (.). Hereafter, square matrices in 
the replica space are denoted by a double overbar. The identity in the replica space is J" 7 = S ay , and we furthermore 
define the matrix U ai — 1 for all a, 7, so that 



hi - (3(v/V)M 2 U. 



(3.11) 



The second-order term in the weak-disorder expansion of (Z^ t ) is provided by the gaussian integral ([T^) only, 
neglecting the corrections 0(e 3 ) and higher. As in centered scalar gaussian integrals, the term proportional to e 3 
vanishes upon integration, and the next non-zero correction in fact is an 0(e 4 ). The aim of the paper being to explain 
the functional method, we shall ignore such higher-order terms hereafter but they could as well be computed by 
pursuing Laplace's expansion scheme. 

The details of the integration over v with measure T> ut v are left to Appendix Taking into account that e is 
expressed in terms of the gradient of v involves a change from v(x) to its Fourier components v(k). One finds 



(^}oc{e- r W ) ^n dct [Q £ . (k.M.fc) .Q 

I 2 [.Lit], rep L v 



-1/2 



(3.12) 



where k = k/fc and Q £ = I — kk. 

We now calculate the determinant. Defining 

Mi = Q fc .(k.M 1 .k).Q £ , 
M 2 = Q £ .(k.M 2 .k).Q £ , 

we obtain, with the help of formulae (C2) in Appendix ^|: 



det 

2[_Lk],rep 



k .Qt 



log det (Mi) 

2LLkl 



/?-^tr(M 2 .M~ 



0(r 2 



(3.13a) 
(3.13b) 



(3.14) 
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Carrying out as in ( 2.1 8| ) the limit r — > in (3.12) it follows that 

1 ,. , . , 1 v 



+ o(log(/?)//3) + 0(1/0). (3.15) 

The sum over k ^ . . . (divided by V) can now be replaced by a <i-dimensional integral [divided by (27r) d ] over all the 
Fourier modes k in the continuum limit, since the integrand is not singular at k = 0. The finite size of the elementary 
cells of volume v is accounted for by an upper cut-off A in this integral. It is chosen such that 

•Lw- 1 (3 - 16) 

Hence A -» oo when v -> 0. We denote by S d = 2n d / 2 /T(d/2) the area of the sphere of unit radius in R d , and by 
dOj^ the integration measu re over the solid angl e in the direction k (dk — dkdft^). 
Letting f3 — » oo in |Tl5| ) we deduce, by ( pl2| ): 



$ut(D) - <0 X (D)) - ~ / ^ tr (m 2 .M x 12 



2] S d 2 

This result is the second-order perturbation expansion valid for any potential </> x (d) with the constraint tr 2 d = 



(3.17) 



2. Application to the Norton law 



We evaluate the previous result ( $-17 ) for the Norton potential (|3.lD. We write 



.)+■ 



Defining D = D'/ v /tr 2 (D' 2 ), so that tr 2 (D.D) = tr 4 (DD) = 1, (|3/]) and (gj) entail: 

Mi = mil + to 2 DD, 
M 2 = to 3 DD, 

where 



d ^""" e q ' 
m 2 = (to — 1)toi, 

Setting 

£) t = |Q fc .D.k| 2 = k.D 2 .k- (k.D.fc) 2 , 
and introducing the unit vector u = Qg..D.k/|Qj..D.k|, we obtain 



Mr = f (Q k 



i) + (m 2 ^ £ 



TOl 



M, 



(3.18) 

(3.19a) 
(3.19b) 

(3.20a) 
(3.20b) 
(3.20c) 

(3.21) 



(3.22a) 
(3.22b) 



k r 2 

TO 3 D fc uu. 

Since Qg. — uu and uu, are a pair of complementary orthogonal projectors in_L(2[_L k]) (where Q £ plays the role of 
the identity), the inverse of Mi, its product with M 2 and the final trace in (3.17) are readily obtained, resulting in: 



$ut(D) 



c(D)> 



Sd toi + 2to 2 Z? £ 
Now replacing mi, m 2 and m 3 by their values ( 3.20| ) finally yields 



$ ut ( D ) = gmcff(D) ur +1 



TO + 1 



«) / (m + l)£> £ 3 3 ' 



flmcff(D) = 

This result reduces to that of Suquet and Ponte-Castaneda (1993) for a space dimension d = 3. 



(3.23) 

(3.24a) 
(3.24b) 
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B. The viscoplastic potential 



We undertake a similar cal culat ion for the viscoplastic potential. It can be ch ecked (cf. Appendix E) that the result- 
ing perturbative expression ( 3.34 ) is linked by a Legendre transform to ( 3.24 ) , the one for the dissipation potential, 
though different boundary conditions are considered. This is a consequence of our working in the thermodynamic 
limit where the volume of the system V — > oo (in the calculations, this shows up essentially in that the Fourier sums 
are carried out with k > 0, rather than with k > (2tt)/L, where L would be the typical system size). As noticed 
by Hill and Mandel in elasticity, the homogeneized stiffness tensor obtained with uniform strains on dft, and the 
homogeneized compliances tensor obtained with uniform stresses on <9fi are inverse of one another up to terms of 
order 0((£/L) 3 ) in three dimensions, where £ is the typical size of the heterogeneities |2(|. They therefore become 
exact inverses in the thermodynamic limit. The same effect is at play here. The viscoplastic potential can thus 
directly be obtained from the dissipation potential as done in Appendix E. However we prefer to compute it from 
scratch because differences with the previous section show up, and it may be useful to see what these are for future 
purposes: indeed, as soon as one leaves exact perturbative expansions and deals with self-consistent estimates, one 
approximation scheme worked out separately for dual potentials quite often generates two estimates which are not 
duals of one another (a notable exception is the self-dual EMT theory). Treating in a symmetrical way the two 
potentials is then a more sensible thing to do than working with only one. 



1. General framework 



We now have to compute 



(3.25) 



The perturbative expansion of the effective potential is carried out around the constant solution cr = S, writing 
a = S + s. The expansion of i/j x (E + s) has the same form as (3.6). We now denote by W and W the first and second 
derivatives of tp x (cr) evaluated at cr — X, respectively. Using s as integration variables, the expansion of (Z^ s ) is: 



( z us) cc (e-W**™)* I ]Jv s s a [l + 0(s 3 )] 



x exp 



(3.26) 



where V s s — D s s<5(V.s)5(s) and MP 7 is defined as in Eqs. (3.9) and (3.11). Once again, the details of the integration 
over s with measure T> s s are left to Appendix The implementation of the constraint V.s = involves a change 
from s(x) to its Fourier components s(k). The result is: 



(Z r na ) cx ( e - r ^WS)\£ TT J det ( M ) det 



4s', rep 2[J_k],rcp 



4s ' rop .k).Q t 



-1/2 



(3.27) 



The determinants are once again straightforwardly obtained, and the limits r — * and (3 — > oo carried out as in 
the previous section. Setting now 



Mi =Q £ .(k.M- 1 -'.k).Q fc , 

M 2 = Q k .(k.M7 l4s ' :M 2 :Mj; lls '.k).Q fc , 

the second-order perturbation expansion for the viscoplastic potential finally reads 

1 f dfl 



* US (S) = (V> x (£)) 



S d 



tr(f 



tr(M 2 .M 1 2[±k] 



(3.28a) 
(3.28b) 



(3.29) 
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2. Application to the Norton law 



Writing w„ = (u n ) + Suj n , and E = E'/ ^/tr 2 (E' 2 ), we have 



11 = mj + m 2 SE = tox (J - EE) + (mi + m 2 )EE 

1 2 = to 3 EE, 



where 



(3.30a) 
(3.30b) 



TO l 



J 1 ( W ")^eq S 



777-2 = (n — 1)7771 



7773 



2\ v>2n 
eq • 



(3.31a) 
(3.31b) 
(3.31c) 



Noting that JJ — EE and EE are complementary orthogonal projectors in L(4s'), we compute M. 1 l4s ' , M. 1 li3 ' : M 2 
and M 1 1<ls ' :M 2 :M 1 l4s '. Now introducing the unit vector u = Qg..E.k/|Qg..E.k|, we express Mi and M 2 in terms of 
Qr — uu and uu and compute the rest. With 



S c = l-2|Q t .E.k| 2 = 1-2 



A 2 „- 



k.E .k- (k.E.k)' 



we arrive at 



^(E ) = (^ ( E))-i/^-^. 

2 J S d toi + to 2 E £ 



(3.32) 



(3.33) 



Finally replacing mi, to 2 , 7773 with their values (3.31) yields the desired expression for the effective Norton viscoplastic 
potential: 



* US(S) " 77+1 ^ ' 



W ne ff(E) = (w n ) 



5 d l + (n-l)E 



(3.34a) 
(3.34b) 



IV. OTHER TYPES OF DISORDER 



We briefly discuss here how to cope with non site-d isordered models, with possibly anisotropic correlations. It 
suffices to write the average of the exponential in ( p. 20 ) simply as 



-/3E;=i0x(d-r)\ = 



exp 



log/e-^W^) 



(4.1) 



and to expand the logarithm as in Sec. Ill A 1, only keeping the quadratic terms in the outer exponential. Equ. ( |3.8| ) 
is unchanged, save for a slightly different (but equivalent) form for the first averaged prefactor, and for a modification 
in the quadratic term in the main exponential. This term becomes 



2 V 



£ J dxdye Q (x):M^(x|y):e Q (y), 



(4.2) 



where 



rra~f 



(x|y) = (W")<W(x - y) - | [(W'(x)W'(y)) - (W'(x)> (W'(y))] . 



(4.3) 



This requires modifications in the calculations, but does not change their principle. One merely has to deal with 
summations over another pair of indices, namely the coordinates x and y. Now the determinants, e.g. in (3.12), 
become functional determinants and must simultaneously be evaluated as determinants of operators in the Fourier 
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space, in addition to being operators on finite-dimensional vector spaces. However, this is not a serious problem in 
the perturbative calculation considered here since only traces are required in the final second-order results. 

In the case of translation-invariant disorder, the kernel M" 7 (x|y) is diagonal in the Fourier representation and 
( 3.12 ) still holds, provided that one writes M(k) (for its Fourier transform with respect to x — y) instead of M. The 
two-point correlation function g(x) is defined by: 



(Suj n (x)Suj n (y)) = (Sul) g(x - y), 



with 



S(x = 0) 



dk 



.9(k) = 1. 



The final results are then only modified as 

? 2 ) 

m / 



cff(D) = 



1 - 



dk (m + l)s(k)D fi 



W„ e ff(S) 



1 



1 (Sv 2 n) 



(27r) d 1 + 2(m - 1)D £ 

dk (n + l)g(k)S fc 
(27r) rf 1 + (n - 1)E£ " 



0(W/(^m) 3 ) 



0«^)/( w „) 3 ) 



(4.4a) 
(4.4b) 

(4.5) 

(4.6a) 
(4.6b) 



The Legendre duality property once again applies, whatever the nature of the correlations, and the result is once 
again independent of the boundary conditions. The previous expressions for isotropic site-disorder are recovered using 
g(x) = vS(x). Ellipsoidal symmetry corresponds to a correlation function g(x) = g(|Z.x|), where Z is a symmetric 
constant tensor (Willis 1977; Ponte-Castaheda and Suquet, 1995). 



V. DISCUSSION 



In this section, we address the range of applicability of the functional method. Firstly, a major concern of studies on 
heterogeneous media is about bounds. In prin ciple, the a bove functional formulation should be apt to deliver bounds, 
at least as far as the initial expressions (2.12) and ( 2.15| ) are concerned. One might think, e.g. to use the convexity 
properties of the exponential in order to obtain inequalities on approximations to the potentials. However, to the 
knowledge of these authors, the replica method is an inadequate tool to use for subsequent calculations, if these are 
to preserve inequalities. Indeed, a convexity equality of the type 



(*-v) 



> e 



'(i:; = i v=) 



(5.1) 



where [i is some normalized functional measure and V is some potential, valid as long as the potential is replicated 
r times, r > 1, may not survive the limit r — > 0. In the zero-replica limit indeed, quantities always positive when 
r > 1 may become negative, minima may transform into maxima (Parisi, 1984), and it is most often hard to conclude 
on inequalities. The replica method has been invented, and employed, to find estimates (which can be of very good 
quality) to free energies when disorder is present, not bounds. 

Second, one must not confuse the initial functional formulation with the replica method, which we use as a tool in 
order to work out averages. These are distinct matters. Actually, there are two starting points for the calculations: 
either we transfer the statistical averages of the logarithm inside the functional integral by means of the replica 
method, and carry out approximations on the effective potential 

-Ilog(e-^^)), (5.2) 



cf. (4.1). This is what we did in the paper. Either we carry out an approximation on the non-replicated initial 
functional integral ( 2.12 ), then take the logarithm, and average afterwards only. Save for perturbative expansions, 
the results will be different. 

Third, all types of calculations feasible with classical m etho ds (es pecia lly the perturbative ones) should translate 
into simple approximation schemes to the basic equations (2.12) and (2.15). But we believe that the compact starting 
point of the functional formulation might enable one to obtain non-perturbative results more easily, and of different 
nature, than by classical means. There resides its main interest. 

Finally, the trick consisting in writing the minimum of some functional as the limit of a particular functional integral 
could be applied elsewhere, e.g., to compute variational expressions such as those presented in Ref. p3| (computing 
a maximum can of course be done by identical means). 
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VI. CONCLUSION 



A field-theoretic method in order to compute the effective properties of highly non-linear composites has been 
introduced. The minimization problem was shown to be equivalent to the computation of a functional integral. 
We believe the main advantage of this formulation is that it is valid for all types of local potentials, and that it 
can be used as a convenient starting point for approximations. As a first application, a weak-disorder second-order 
perturbative calculation of this functional integral was presented. The second-order results of Suquet and Ponte- 
Castaneda (1993) for the effective plastic dissipation potential (with uniform traction boundary conditions) were 
recovered. In addition, the second-order correction to the viscoplastic potential was also obtained, and these results 
were extended to correlated disorder. With this method, minimizing over a vector velocity field, or a tensor stress 
field, is found to be equally feasible. 

A natural extension of this work is to find a well-behaved non-perturbative approximation. Self-consistent calcula- 
tions are currently under study. 
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APPENDIX A: FOURTH-RANK TENSORS IN MECHANICS 



Before briefly reminding the reader about the main properties of fourth-rank tensors in mechanics, we need to set 
our notations about the various algebra and sub-algebras of linear operators that come into play. 

We designate by L(2) the algebra of linear operators (second-rank d x d square matrices) M : R d — > R d . Next, 
fourth-rank tensors appear in mechanics as representations of linear operators M : L(2) — > L(2). They can therefore 
be considered as square matrices, the indices of which are pairs of indices, and their components are denoted by Mij t ki, 
such that i, j, k, I = 1, . . . , d. We denote the algebra they belong to by 1/(4). The product in L(A) is 

(M:N)y >fc | = M l0 . mn N nmM . (Al) 

Hence applying M to a second-rank tensor N , we have 

(M:N) tf =M ijM N lk . (A2) 



The identity I in 1,(4) is 

and the trace is defined (and denoted) by 



hj,ki = SuSjk, (A3) 



tr 4 M = M ijtji . (A4) 



Hence tr 4 I = d 2 . 

The algebra L(2) can be split in the direct sum of that of symmetric traceless matrices [which we denote by L(2s')], 
diagonal matrices proportional to the identity [L(2d)], and anti-symmetric matrices [L(2a)]. The algebra of symmetric 
matrices is the direct sum L(2s) = L(2s') © L{2d). The algebra L{4) therefore admits three remarkable sub-algebras, 
which we denote by L{As'), L(4d), L{Aa), generated by the mutually orthogonal projectors J, K and I a respectively, 
defined by 

Jij,ki = -^{SuSjk + SikSji) - -Sijdki, (A5a) 
Kij.ki = -jSijSki, (A5b) 
ItoM = \{5ii$ok-kk5ji)- (A5c) 
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They obey I = I a + J + K. Hence, the sub-algebras L(4s'), L(4d), L{Aa) are that of the endomorphisms on L(2s r ), 
L(2d), L(2a) respectively. The operators JJ, K, and 1° are the appropriate identity operators in each sub-algebra. We 
have 

tr 4 J = d(d + l)/2- 1, (A6a) 
tr 4 K = 1, (A6b) 
tr 4 I° = d(d- l)/2. (A6c) 

These numbers count the number of eigenvalues of each operator. In L(4s), the sub-algebra of the endomorphisms 
on L(2s), the identity is P = J + K: 

IliM = ^iiSjk + SikSji), (A7) 

and tr 4 F = d(d + l)/2. Note that L(4s') © L(Ad) C L(4s), but the two sets are not equal. 

Finally note that the inverse (resp. determinant) in L(4) of an operator M S L(As') © L(Ad) © L(4a) is the sum 
(resp. product) of the inverses (reps, determinants) of its constituents in their respective sub-algebra. 



APPENDIX B: GAUSSIAN INTEGRALS 



1. Gaussian integrals over vector fields 



We start from the well-known d-dimensional gaussian integral over a vector field A (Kleinert, 1995). If M is a d x d 
positive-definite symmetric matrix, and b any vector, then 



i/A cxp j — ^A.M.A + b.A 



(27T) 



det 2 (M) 



1/2 



exp Qb.M~ l2 .b^ , 



(Bl) 



where the integration measure is d\ = \\ i d\i . 

This formula is meaningless if M is not invertible in L{2). However, let k be a unit vector, b now be a real vector, 
and M = M± + e 2 kk, where M_l is invertible only in L(2[J_ k]), the set of the endomorphisms which operate in the 



subspace of vectors orthogonal to k. Then 



MJ 12 ^' +e- 2 kk, det 2 (M) =e 2 det 2[±£] ( 



and 



d\ exp [ --A.M. A + ib.A 



(27T) 



e 2 det 2fJ _ fcl (M^) 



1/2 



exp 



.b-— (b.k) 2 



Letting e — > and using the representation of the Dirac distribution 



5(x) = lim 



1 



|(x/e) 2 



we obtain 



/ 



<l\ ck)> j -^A.IVL.A + ib.A 



o V^ 2 



(2 7 r)' i + 1 



1/2 



_det 2[xfc] (M ± ) 

We define the integration measure over the vectors orthogonal to k to be 

dX = S(\.k)d\. 

[±k] 



exp (^-^b.M^^.b^ 5(b.k) 



(B2) 



(B3) 



(B4) 



(B5) 



Multiplying both sides of (B4) by exp(ib.a), where a is a real vector, and integrating over b yields 
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/ 



db exp ( — ib 

[Xfc] 



(27T) 



-(d+1) 



n 1/2 



= (27T) d 



r -,CM _ 2[J - kl 

ut!u 2[Xk]V m ± 

(27T)-( d+1 ) 



b + ih.a. 

dXdh exp 
-,1/2 



— A.Mj_.A 

2 



det 2[±£] (r 
(27r) d " 1 







1 1/2 



det 2u£l (M 



-lar. 



exp 



dX exp ( --A.Mj..A 



— a.Mj_.a 
2 



zb.(A + a) 
S(a + A 



where <5(x) = J(xi). Note that we can exchange M 



and Mj_ in ( 



Moreover, 



(B6) 



can be written 



Mj^ = Qk-M.Qf., where Q £ is the projector Q £ = I — kk, for some M 6 L(2). Besides, the b are now orthogonal to k 
(because of the integration measure) so that b.Q £ .M.Q £ .b = b.M.b. Finally, the result can be analytically continued 
to complex values of a. Changing the names of the variables, we therefore arrive at the formula 

1/2 



/ 



dX exp 

[Xt] 



1 



A + b.A 



(2vr 



,d— 1 



det 2[J . £] (Q fc .M.Q £ ) 



exp 



ib.(Q £ .M.Q t )- 1 »ixs ] .b 



(B7) 



2. Gaussian integrals over tensors of rank 2 



The generic gaussian integral over all matrices s £ L(2) is 



ds exp 



1 

2 S: 



b:s 



(27T) 



det 4 (M) 



1/2 



exp 



(B8) 



where ds = Yin dstj is the appropriate measure, M is a symmetric matrix of £(4), i.e. is such that My^t = Mki,i 



This formula is a direct consequence of (Bl) via a mapping of s onto a column vector in R rf . Once again, the integral 
is meaningless in general if M is not invertiblc in L(4). We apply the same method as in Sec. B 1. 



a. Gaussian integrals over symmetric matrices 



Let b £ L{2) be real. We consider M = M s + (e 2 /2)P, where M s £ L(4s) and is invertible in L(As). Then 
I" 1 " = Mj l4s + (2/e 2 )F and det 4 (M) = det 4s (M)(e 2 /2) rf ( d - 1 )/ 2 . Hence from @, 

1 



ds exp ( — -s:M:s + ib:s 



det 4s (M s ) 



1/2 



( e 2/ 2 )d(d-l)/4 



'' "I7 l4s :b- 4rb:I a :b 

e 2 



But b:I a : b = (1/2) Ei<, (&ij " b ji) 2 - Using @ and going to the limit e -> 0, we find 



e?s exp ( — s : M s : s + ib : s 



1/2 



exp(--b:M; 1 -:bm^-^) 



i<j 



(2ir) d2 (4ir) d ( d -V/ 2 
det 4s (M s ) 

We define the measure over the symmetric tensors of rank two as: 

dsS = 2 ^- 1 )/4[-Q S ( Sij - Sji )]ds. 

i,j<i 



(B9) 



(BIO) 



(Bll) 
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Multiplying both sides of ( B10| ) by exp(ib:a), and integrating with respect to b with db yields 

1. 



d s b exp ( -^b:M7 l4s :b + ia:b 



(27T) 



-d 2 -d(d-l)/2 



det 4s (M^ l4s ) 



1/2 



1 



= (27T) d 



(2 7r )- d2 - d ( d - 1 )/ 2 



dbds exp I — -s:M s :s + i(s + a) : b 
1 1/2 



det 4s (M7 l4 

( 27r )d(rf+l)/2 1 1/2 



ds exp ( --s:M s :s ) <5(s + a) 



det 4s (ML7 i4s ) 



exp ( --a:iw s :a 



(B12) 



where <5(x) = J}- <5(xjj). Whence the generic result for gaussian integrals over symmetric matrices, with M G L{A): 



il^s exp j — -s:M:s + b:s 



(• 27r )d(d+l)/2 



1/2 



. .xj> ( ib:(F:M:I s )- l4s :b 



_det 4s (I s :M:P) 
6. Gaussian integrals over traceless symmetric matrices 



(B13) 



is not invertible in L(4s), but only in L{As') for instance, the procedure can be repeated: let us assume 

;.p)-i4, = M7, l4s ' + K/e 2 , 



If F:R 
that M = 

det 4s (I s :M:I S ) = det 4s /(M s /)e 2 and letting e -> 0, one findsT 



e 2 K, where M s / G L(As'). Starting from (B13) with b — > ib, using (I 



d s s exp ( — -s:M S ' :s + ib:s 



d(27r )d(d+l)/2+l 



1/2 



deW(M s ,) 

We define the integration measure over symmetric traceless tensors to be 

d B >s = d 1/2 8(tr 2 s) d e s. 



exp ( -ib:M s > 3 ' :b ) S(tr 2 b). 



As one easily checks, we have, with (Bll) 

J d s x exp( lS :x) = 2 d ( d - 1 )/ 4 (2^) d ( d+1 )/ 2 [[] ( 5( Sli )][J]5(^ + Sji )}. 



i<j 



Multiplying (B14) by exp(ia:b), and integrating over b with measure d s b, we deduce, using 



for gaussian integrals over symmetric traceless tensors: 



/ 



d s is exp [ — -s:M:s + b:s 





"(2 7r ) £i ( d + 1 )/ 2 - 1 " 


1/2 


)- 


_deW(J:M:J)_ 


exp 



-i**':b 



(B14) 



(B15) 



(B16) 



the result 



(B17) 



c. Other gaussian integrals 

For the sake of completeness, we finally give without demonstration the results for gaussian integrals over antisym- 
metric tensors, with measure 



d a s = 2 d ^/ 4 [Y[S( Sll )][Y[5(s l3 + Sji )]ds, 



and over diagonal tensors proportional to the identity, of the type s\/d, with measure 

d d s = d- 1/2 ds J ds6(s-s\/d). 



(B18) 



(B19) 
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One finds 



d a s exp | --s: 



d d s exp [ --s: 



:s + b:s 



:s + b:s 



( 27r )d(d-l)/2 

det 4 a(I a :M:I a ) 
2tt 



1/2 



exp 



_det 4(i (K:M:K) 
The last identity is trivial (scalar gaussian integral). 



1/2 



exp 



-b:(F:I 



-b:(K:M:]K)" l4d :b 



(B20) 
(B21) 



APPENDIX C: INVERSE AND DETERMINANT IN REPLICA SPACE 

In replica space, A and B being replica-independent operators pertaining to some algebra A, the inverse and 
determinant of a matrix of the type 



read 



M 



M = AI + BU, 



A~ lA I -{A + rB)- lA BA~ lA U, 



det (M) = det(A) r ~ 1 det(A + rB) = 1 + r log det (A) + tr^A- 1 ^) + 0(r 2 ) 

A,rep A A L A A 

The expansion derives from the operator identity (Kleinert, 1995) log det = trlog. 



(CI) 

(C2a) 
(C2b) 



APPENDIX D: GAUSSIAN INTEGRALS OVER THE VELOCITY AND STRESS FIELDS 

1. Integration over the velocity field 



We detail here the steps leading to (3.12). 

The Dirac distribution implementing the constraint V.v = in ([T^) is exponentiated first: we introduce a scalar 
field A and write: 



6 (V.v) oc J VX exp i J dx A(x)V.v(x) 



(Dl) 



Then we Fourier transform the integrand. By definition, for any function / (which will be A or s), the Fourier 
transform is: 



/(k) = J dx e lk x /(x), /(x) = /ct e- ik '*/(k). 
Since /(x) is real, we write for k ^ 0: 

/(k) = i=[/W(k) +i /( 2 )(k)], 
/(-k) = i=[/W(k)- i /( 2 )(k)], 



(D2) 

(D3a) 
(D3b) 



where and are real functions of the wavevector; also, we set f(k = 0) = /W(0). This decomposition is 
carried out for all A" (k) and v Q (k) . 

Introducing the set K. + = {k/fej > 0} U {k/ki — 0, k 2 > 0} U . . . U {k/ki = 0, . . . , fc d _i = 0, k d > 0}, a functional 
measure T>f over /(x) thus becomes in the Fourier representation 



Vf = n#(x) d/W(0) [] Wf (1) (k)d/ (2) (k)]- 



(D4) 
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Note that /W and are denned only for k £ /C + . By Parseval's identity, we have 



dxA(x)V.v(x) 



(2tt)° 



A< 2 > (k)v« (k).k - A« (k)v( 2 ) (k) .k 



(D5) 



As usual in functional methods, the latter integral is to be understood as a Riemann discrete sum by applying the 
correspondence: 



dk 



(27T) 



d y 



(D6) 



I n p articular, this allows one to separate the contribution of the mode k = 0. The argument of the exponential in 
( |3.8D is likewise transformed into a sum of Fourier modes via 

J rfxe Q (x):M Q7 :e 7 (x) 

= I e a W (o) : M q t : e^ 1 ) (0) + - ^ [e"' 1 ' (k) : MP 7 : e 7 « (k) + e a ^ (k) : M q t : e^ 2 > (k) 

kGK:+ 

= i J2 [v a(1) (k).(k.M a ' 1 '.k).v 7 W(k)+v a ( 2 )(k).(k.M a! 'T.k).v^( 2 )(k) 



(D7) 



In the last equality, we have used the symmetry of M Q7 with respect to its tensor indices. Note that e^ 1 ' = — [kv^ 2 ) + 
v( 2 )k]/2 and e< 2 ) = [kv« + v«k]/2. 

Because of the constraints S(e a ), the linear terms e Q disappear from the exponential. There is no dependence of 
the integrand with respect to AW(0) nor to vW(0), so that integration over these variables only yields harmless 
/3-independent infinite multiplicative factors which we drop out, since they do not contribute to the final result in the 
limit P — > oo. The reason is the same as that invoked in conjunction with Equ. (2.24). Such physically irrelevant 
(because multiplicative) infinities are often encountered when dealing with functional integrals. They are related to 
(here) unimportant normalization questions. 

We therefore see that for each k ^ 0, mutually complex conjug ate partial integrals A(k) = J d\ a ^(k) dv Q W(k) 
[. . . ] and A* (k) = / d\ a(1 > (k)dv a W (k) [. . . ] show up in pairs, and can be evaluated independently from one another. 
Renaming, in A(k), the integration variables into k-independent d\ a and d~v a , we arrive at 



oc <e- r "**«< D >> - J] \ A( ^)\' 



where 



exp 



A(k) = I I H dA " rfv 

\ a / 

J[dv a 5(v a .k) 

a 

dct [Q £ . (k.M.k) .Q fc 

2[_Lk],rep L \ ' 



- \4* E v " • ( k - MQ7 - k ) - v7 + h E AQ ( vQ - k ) 



2V 



exp 



c*7 
-1/2 



^v a .(k.M a ' lf .k).v' ) 



(D8) 

(D9a) 

(D9b) 
(D9c) 



and Q £ = I — kk with k = k/k. In the step from ( D9b ) to ( D9c ) we used formula (B7), triviall y gen eralized to the 
replica space, and discarded, among other irrelevant factors, a power of the modulus k. Equ. (3.12) finally follows 
from the fact that A(k) = A(-k). 



2. Integration over the stress field 



We now detail the steps leading to ( 3.27 ), which differ little from above. Because of our considering the Norton 
viscoplastic potential (which does not depend on <r m ) for explicit applications, the matrix Mi given in (3.30) is not 

invertible in L(4s), but only in L(4s'). As a consequence [cf. (C2a)], M is not invertible in L(4s,rep), but only 



17 



in L (4s', rep). This leads to peculiarities that have to be taken into account but does not change the principle of 



the calculation. The equivalent of formula ( 3.27 ) for M invertible in L (4s, rep) is provided hereafter for the sake of 
completeness. 

The Dirac distribution implementing the constraint V.s = in the measure T> s s is exponentiated by the introduction 
of a vector field A as: 



(D10) 



<5(V.s) oc J V\ expfi |dxA(x).[V.s(x) 

= ^Aexpfii {s (1) (k):[kAt 2 Hk)]-s( 2 >(k):[kAW(k)]} 



ke/C-| 



A reasoning paralleling that in Sec. D 1 yields 



VEK+ 



where now 



A(k) = J \ \\d\ a d s s a J exp 

-1/2 



1JL 

2V 2 



s" :M Q1 :s 7 + ^ ^ s a : (k A Q ) 



oc det 

4s' ,rep 



cm / ^ ay 



oc det (Mfp 1 / 2 det 

4s', rep 2[J_k],rep 



Q t . (k.M l4 *'- rBP .k] .Q t 



k. fM~ 1 *«'.«*') a ' r .k 



-1/2 



(Dll) 

(D12a) 
(D12b) 
(D12c) 



and Qj- = I — kk, and k = k/fc. In the first step from (D12a) to (D12b), we used the formula (B14) generalized to the 



replica space. The resulting factor <5(tr2 kA) oc <5(k.A) was absorbed in the measure d\ a ^ ± ^ defined in (B5) 



Next, 



we appealed to (B7), also extended to replica space. Equ. ( 3.27 ) follows. Had M been invertible in L(4s,rep), one 
would have found 



Z r ) oc e~ rp v 



r/3fV-(E)\ 



IT \ det (M) det 

I 4s, rep 2, rep 



- l4s.rep -v 

k.M ' .k 



-1/2 



(D13) 



instead. 



APPENDIX E: THE LEGENDRE TRANSFORM 



Though the plastic dissipation potential (3.24) and the viscoplastic potential (3.34) have been obtained for different 
boundary conditions, we show here that they are linked by the Legendre transform ( [2.8] ). 

Let us deduce (3.24) from ( 3.34| ), for instance, assuming that ( |2.8| ) holds between both. We write $(D) = $o(D) + 
<5$(D), *(E) = * (£) + <$ *(£), where $ (D) = (<MD)), *o(S) = (V'x(S)) are the leading terms in the perturbative 
expansions, and <5$(D) and ^(E) are the corrective terms. It is straightforward to check that $o(D) and ^o(S) are 
Legendre duals: $ (D) + *o(S) = E : D. Hence 

$(D) = E:D- * (E) -«y*(S) = $ (D) - £*(E), (El) 

so that (5$(D) = -<5*(E). Moreover, since E = <9$(D)/<9D = d$ {D)/dD + <9<S$(D)/<9D, we obtain to second order 

$(D) ~ $ (D) -5*(a$ (D)/aD). (E2) 



The Legendre transform of V&o(S) = (w„)S"+ 1 /(n + 1) reads, with 



l/n: 



$o(D) = (wn)-™ U2 + V(m + 1) 



tr(D) = 0. 

2 



(E3) 
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In addition, 



d- 1 



-'cq — x^n 



(E4) 



These expressions for S' and S eq are correct to first order only. However, symmetry considerations show that S' is 
always proportional to D. Thus X = D, whence from the definitions ( 3.2l| ) and ( 3.32p , 



The weak-disorder expansion of ^(D) has to be expressed in terms of m = 
we have 



-l/m 



im + i (sel 



2 m 2 



(E5) 

To second order, 
(E6) 



Gathering these results into the perturbative expansion (E2) of $(D) finally yields (3.24) back, as announced. 
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